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We study mesoscopic resonant tunneling as well as multichannel Kondo problems by mapping 
them to a first-quantized quantum mechanical model of a particle moving in a multi-dimensional 
periodic potential with Ohmic dissipation. From a renormalization group analysis, we obtain phase 
diagrams of the quantum Brownian motion model with various lattice symmetries. For a symmorphic 
lattice, there are two phases at T = 0: a localized phase in which the particle is trapped in a potential 
minimum, and a free phase in which the particle is unaffected by the periodic potential. For a non- 
symmorphic lattice, however, there may be an additional intermediate phase in which the particle 
is neither localized nor completely free. The fixed point governing the intermediate phase is shown 
to be identical to the well-known multichannel Kondo fixed point in the Toulouse limit as well as 
the resonance fixed point of a quantum dot model and a double-barrier Luttinger liquid model. The 
mapping allows us to compute the fixed-point mobility /i* of the quantum Brownian motion model 
exactly, using known conformal-field-theory results of the Kondo problem. From the mobility, we 
find that the peak value of the conductance resonance of a spin- 1/2 quantum dot problem is given 
by e 2 /2h. The scaling form of the resonance line shape is predicted. 
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I. INTRODUCTION 



Electronic transport in mesoscopic systems has drawn much attention in the last decade due to recent development 
of fabrication techniques of ultra small structures of high quality. In the nanometer scale systems, electron-electron 
interactions play more important roles, as the number of effective spatial dimensions d is lowered,-. Most remarkable 
among many examples are fractional quantum Hall effeet and Wigner crystallization at d — 2, El Luttinger liquid 
behavior at d = Ipu and Coulomb blockade at d = O.B On the theoretical side, the use of traditional perturba- 
tive approaches were proven less useful due to the strong correlations, while nonperturbative techniques such as 
renormalization group, conformal field theory, and the Bcthe ansatz have seen some success. 

Studying resonant tunneling in a quantum dot coupled to leads via two identical quantum point contacts (QPC), 
Furusaki and MatveevQ showed that at T = CLthc on-resonance conductance of a spin- 1/2 system takes a universal 
value less than the perfect conductance e 2 /hx\ They argued that the on-resonance properties are controlled by the 
highly nontrivial non-Fermi-liquid fixed point of a four-channel Kondo problem. The deep connection between the 
problem of resonant tuancling in a point-contact-coupled quantum dot structure and the -Kondo problem has been 
known for some time.u'u A recent experimented has confirmed the theoretical predictions" of the capacitance line 
shape which was based,on this connection. 

Later, Yi and KaneE3 showed that there is a general and direct one-to-one mapping between the resonant tunneling 
problem and the Kondo problem in a special limit. They argued that this limit is a multichannel generalization of 
the well-known Toulouse limit .113 Using the conformal- field-theory results on the Kondo fixed point, they computed 
the universal conductance on resonance exactly. 

They have also found that the quantum-field-theoretical formalism of both problems can be reduced to a first- 
quantizcd-.auantum mechanical formalism of a fictitious particle that moves in a periodic potential in a dissipative 
medium.EjhJ Technically, this single-particle quantum Brownian motion (QBM) model is obtained by integrating 
out all degrees of freedom in the original quantum field theory except only a small finite number of them. These 
remaining degrees of freedom become the "spatial coordinates" of the Brownian particle and the integrated ones serve 
as the source of effective dissipation. 

The symmetry of the periodic potential in the QBM model has been found to be crucial in identifying the zero- 
temperature phase diagram, while a particular phase is determined by the potential period and the dissipation 
strength. In particular, the on-resonance quantum dot problem were related to potentials with non-symmorphic 
lattice symmetries. A similar QBM model has been used |-i«|-*elatively recent studies in the context of resonance 
tunneling through an impurity state in a Luttinger liquidJIU'ta which apply to such experimental situations as in 
quantum wiresoO and quantum Hall edge systems.B More applications of the QBM descriptions may also be found in 
connection with tunneling between quantum Hall edge states, both in fractionalcJ and integerEH quantum Hall regin 
A similar approach has also been used to study resonant point-contact tunneling between Luttinger-liquid leads, 
but the QBM model in that work was formulated with a 7r-flux through each plaquette. 

In this paper, we extend the investigations in Ref. ^ and present a more general description of the QBM model as 
applied to (i) resonant tunneling in a quantum dot coupled to an arbitrary number of leads, (ii) multichannel Kondo 
problem, and (iii) resonant tunneling of Luttinger liquid through a double-barrier. As will be shown below, the QBM 
model possesses very rich phase diagrams. There are three zero-temperature phases: (i) the "free" phase in which the 
Brownian particle completely ignores the periodic potential, (ii) the "localized" phase in which the particle is trapped 
in a minimum of the periodic potential, (iii) and the most exotic of them all, the "intermediate" phase in which the 
particle is neither localized nor completely free. Each phase is characterized by a "mobility" fj, — the ratio of the 
average velocity of the Brownian particle to the driving force in the linear response regime — which is universal in 
the sense that it is independent of the potential strength or other microscopic details. The mobility takes a perfect 
value in the free phase (/x = 1), vanishes in the localized phase (/i = 0), and takes a finite value between the two in 
the intermediate phase (0 < /i < 1). 

This paper is organized as follows. Sec. || is devoted to a general and qualitative overview of the QBM description 
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of resonant tunneling and the multichannel Kondo problem. In Sec. Ill, we show how the one-dimensional (ID) 



QBM model is derived from the two-lead quantum dot model of spinless electrons and give a brief discussion on the 
universal scaling. We generalize this discussion to the multi-lead quantum dot model and the multichannel Kondo 
model in Sec. IV. Sec. [V]is devoted to the double-barrier resonant tunneling model of a Luttinger liquid. We discuss 
the similarities and differences between this model and the quantum dot model. In Sec. [v| we develop a general 
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formalism of the QBM model and present the dual theories in the weak and strong backscattering limits 
the mobility and introduce the renormalization group flow equation in both limits. In Sec. VII 
theory to specific lattice potentials. The zero-temperature phase diagram is derive d for each lattice 



We define 
we apply the QBM 
It is also shown 



that there is a stable intermediate phase for non-symmorphic lattices. In Sec. VIII , we explicitly show that the QBM 
model is equivalent to the Toulouse limit of the multichannel Kondo model. We also show that the perturbation 
calculations of physically relevant quantities — the scaling dimension of the leading irrelevant operator and the fixed 
point mobility — are consistent with the exact solutions derived from conformal field theory results. In Sec. IX, we 
calculate the on-resonance conductance matrix G a b of the multi-lead quantum dot model in terms of [L. Especially for 
the two-lead spin-1/2 problem, we calculate the conventional source-drain conductance G* on resonance, as well as 
the scaling form of the resonance line shape. Finally, concluding remarks and discussions are given in Sec. [xj Details 
of some lengthy calculations are given in appendices. 



II. QUALITATIVE OVERVIEW 

In this section, we will use a few simple examples to describe qualitatively the mapping between the QBM model, the 
resonant tunneling problem, and the multichannel Kondo problem. Rigorous derivations as well as more complicated 
systems will be discussed in the following sections. 

The mapping between the QBM model and a mesoscopic tunneling problem is most easily understood by considering 
a single QPC formed by gates on a two-dimensional electron gas (2DEG), which is schematically depicted in Fig. 
|l|. Experiments show that the linear conductance G through a_jQPC is quantized in units of e 2 /ft,,EJc3 which is a 
characteristic signature of non-interacting ID electron systems.tJ Indeed, a QPC may be thought of as a wave guide 
of an electron wave function, and each wave mode plays the role of a ID channel—Each fully transmitting ID channel 
contributes one conductance quantum e 2 /h to the total conductance of a QPC.ED 

As explained below, it is possible to see how a QPC is described by a QBM model, without going into too much 
detail. As the first example, we will consider a QPC that fully transmits only one ID channel and we will construct a 
QBM model that reproduces the correct quantized conductance. The current is obviously given by / = —ed t Q, where 
Q is the number of electrons in the source 2DEG lead. The Euclidean action may be written as 

S = Sq + Ssourccj (2-1) 

where So is the unperturbed action, and Source is the source term due to voltage drop V across the QPC. The second 
term is straightforwardly given from electrostatics by 

dTeQ{r)V, (2.2) 
where t = it is the imaginary time. Using the Kubo formula, the linear conductance is written as 

= ^^\J d T^ T ^I(r)I(0)) (2.4) 

e 2 f 

= lim — r / dio' ujlo' (Q(uj)Q(-uj')) . (2.5) 
uj^o h\uj\ J 

If we assume that So has a simple quadratic form and the Hamiltonian does not explicitly depend on time, then it 
must have the following form in order to yield the correct quantized conductance: 

S = IJ du;\Lu\\Q(Lu)\ 2 . (2.6) 

This may be thought of as an action of a QBM model written in terms of the spatial coordinate of the Brownian 
particle, Q. ,— . 

A rigorous derivation of the action using the bosonization methodcS reveals that the above supposition is correct- 
In the boson representation, a generic non-interacting free ID electron system is described by a Euclidean actioncj'EZI 

S B = I drdx [{d T 9) 2 + [v F d x 9f] , (2.7) 
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where vf is the Fermi velocity. The boson field 8(x, r) may be thought of as the local phase of a ID charge density 
wave. Especially, it satisfies dtO = 2nl and d x 9 — 2np, where / is the current and p is the number density of electrons. 
If we define the coordinate system such that x = at the center of the QPC, then 



Q(t) = / dxp(x,t) (2.8) 
Jo 

0(x = 0,t) , nn . 

= — + const. (2.9) 

Integrating out 6(x, r) except for x = 0, one can easily see that the effective action is indeed given by Eq. (|2.6D . 

Note that the non-analytic form of So suggests that the Brownian particle is not free. This is due to the dissipation 
provided by 6{x ^ 0) after they are integrated out. The analogy between the QPC problem and the QBM model is 
completed by mapping the current I to the average "final velocity" of the Brownian particle and the voltage drop 
V to the "uniform external force". The dissipation is indeed necessary to prevent dtQ from increasing indefinitely, 
because we know from the Ohm's law that it is finite. Even though the QBM action may have to be modified for 
more complicated systems, the above mappings of fundamental quantities are universal. 

Let us now include backscattering by pinching off the point contact via gate voltage. It is equivalent to placing a 
backscattering barrier at the center of the QPC that suppresses electron tunneling. If the barrier is strong, tunneling 
is almost suppressed and Q is strongly quantized at integers due to the discrete nature of electron number. On the 
other hand, if the barrier is weak, Q is only weakly quantized. However, as long as the barrier strength is nonzero, 
quantization effect is also nonzero. This effect is most simply taken into account by adding a term periodic in Q to 
the action. Since the period is one, the most generic backscattering action takes the form 



71=1 



OO , 

ar 



t. 



i2 V nQ(r) 



(2.10) 



where v n is a dimensionless backscattering amplitude and r c is a short-time cutoff. In the QBM m odel , this describes 
a periodic potential. The potential gets tilted if the external uniform force term Source in Eq. ( |2.2| ) is also added. 
This is schematically illustrated in Fig. |[ Transferring an electron across the point contact corresponds to moving 
the fictitious Brownian particle from a minimum of the potential to the next minimum. If we define the "mobility" of 
the particle fi as the final velocity (dtQ) divided by the external force eV in the linear response regime, it is obviously 
the same as the electric conductance G up to a constant factor e 2 . 

If the barrier is strong, it is more convenient to think of the QPC system as two almost separate leads that are 
weakly coupled by a small tunneling amplitude t. The corresponding QBM model is described by hopping of the 
Brownian particle between "lattice sites" that lie on deep minima of the periodic potential. In the above example, 
the lattice is a simple ID lattice with lattice constant a = 1. 

Another example in which the ID QBM description is useful is the interacting ID electron system known as 
Luttinger liquid. Again, if we focus on the current past a certain reference point, say x = 0, we may obtain a QBM 
model by integrating out all degrees of freedom except Q, which is again defined as the total number of electrons on 
the positive axis (x > 0). If there is a backscattering center such as a point im pur ity at x = , it gives rise to a periodic 
potential in the QBM model. The resulting total action is analogous to Eqs. (2.6) and ( 2.1 Cl| ) , but the Luttinger-liquid 



correlation mod ifies the coefficient in front of So which characterizes dissipation. Since we can restore the coefficient 



1/2 in Eq. (2.6) by rescaling Q, it may be alternatively said that the Luttinger liquid correlation modifies the period 
of the potential. From this point on, we will adopt the convention in which the dissipation coefficient is fixed at 1/2 
whereas the period a is a variable. 

The ID QBM model was studied over a decade-ago in quite a different context, namely, as a possible description 
of a heavy charged particle moving in a metal.tSlla P.articularly, Fisher and Zwergert3 have used a model of Ohmic 
dissipation proposed earlier by Caldeira and LeggettEJ to show that there are two phases at T = 0: (i) For weak 
dissipation, the particle diffuses freely as if the periodic potential were completely absent, and the mobility fi takes 
its maximum perfect value (free phase), (ii) When the dissipation exceeds a critical value, however, the particle gets 
localized in a minimum of the potential, a nd j i = (localized phases). As we adopt the convention in which the 
dissipation coefficient is fixed at 1/2 [Eq. fl2.6| )1, the period is given by a = 1 at the phase .boundary. Later, the 
above result was rediscovered in the context of single-barrier tunneling in a Luttinger liquid. lilltiJ In the Luttinger 
liquid model, the above two zero-temperature phases correspond to the two different regimes of the electron-electron 
interaction which is typically characterized by the Luttinger liquid correlation parameter g. (i) If the interaction is 
attractive (g > 1), electrons can freely tunnel through any arbitrarily strong barrier (free phase), (ii) For a repulsive 
interaction (g < 1), however, even an arbitrarily weak barrier can completely block transport and the conductance 
vanishes (localized phase). 
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Another interesting example that is mapped to aJD QBM model is resonant tunneling of spinless electrons through a 
quantum dot. In particular, Furusaki and Matvee'vu have recently studied a system of a quantum dot that is connected 
to two leads via QPCs. In their model, when an electron tunnels into the dot, each electron in the dot has to pay 
the Coulomb energy ~ e 2 /C, where C is the capacitance of the dot. If e 2 /C 3> T, Coulomb blockade occurs and 
the tunneling is suppressed unless a nearby gate voltage is tuned in such a way that the Coulomb energy is exactly 
canceled by the chemical potential of the dot, in which case, resonance occurs. Another assumption of the model is 
that the dot is large enough that the distance between the point contacts is much greater than the thermal coherence 
length Lt — hvp/ksT. In this limit, an electron never travels coherently between the two QPCs. Therefore, if we 
define Q\ and Q2 as the number of electrons in the first and the second lead, respectively, they are independent except 
through the charge conservation condition 

n + Qi + Qi = const., (2.11) 

where n is the number of electrons in the dot. If e 2 /C ^> v n /r c ,T, charge fluctuation in the dot is suppressed and 
we may safely integrate out n (or equivalently, —Qi — Qi)- Then, the resulting effective action has the form of a ID 
QBM action in terms of Q = Q% — Q\. The current from one lead to the other is —edtQ- 

Resonance is achieved when two charge states n and n + 1 (n = integer) have the same energy. In that case, 
electrons no longer pay the penalty when tunneling in and out of the dot. Thus, Coulomb blockade is lifted and 
current is allowed to pass through the dot down to T = 0. If T > 0, the conductance is finite even off resonance due 
to thermal fluctuations. However, it decreases as the gate voltage is tuned away from resonance. In Ref. @, Furusaki 
and Matveev have exactly calculated the resonance line shape of the conductance for spinless electrons as a function 
of gate voltage and temperature. One of the most interesting features of their result is that the resonance line follows 
a universal one-parameter scaling function 

G(5V G ,T) = G(J^, (2.12) 

where 5Vg is the deviation of the gate voltage from the on-resonance value. Similar results were .obtained for resonant 
tunneling problems in a Luttinger liquidliS or between fractional quantum Hall edge states I 20 !! 31 ! 

Furusaki and Matveev have also argued that the on-resonance behavior of the spinless quantum dot model is 
governed by the fixed point of the two-channel Kondo problem on the Emery-Kivelson linsx3 This limit may be 
thought-<pf as a generalization of the Toulouse limit in the single-channel Kondo problemlLij and may be solved 
exactly £3 In fact, the on-resonance spinless quantum dot model may be shown to be equivalent to a Kondo problem 
on the Emery-Kivelson line. Specifically, it consists of two ID electron channels and a local magnetic impurity, both 
of which have spin 1/2. Later, we will rigorously prove the equivalence, but it may be rather straightforwardly 
visualized via the following one-to-one mappings: (i) the two leads are mapped to the two Kondo channels; (ii) the 
number of electrons Q\ and Q2 are mapped to the z-components of the total spin Sf and S% in respective Kondo 
channels; (iii) the two degenerate charge states of the dot is mapped to the two spin states of the impurity; (iv) 
the tunneling amplitude t is mapped to the exchange interaction coupling constant J; (v) the resonance driving gate 
voltage deviation SVq is mapped to the local magnetic field TL that applies only to the impurity spin; (vi) the voltage 
drop between the leads eV is analogous to the difference in the external magnetic field 7Y e — Hoc between the channels; 
(vii) electric conductance G is mapped to the "spin conductance" defined as ]hn.H 3 -*Hi 9t(<S| — Sf)/(Hg — Woo)- 

The above Furusaki-Matveev mode is nnitp similar to the spinless Kane-Fisher model of resonant tunneling through 
a double-barrier in a Luttinger liquid.EZEj However, there are several important differences. First of all, the Furusaki- 
Matveev model assumes no interaction in the lead. Consequently, the quasiparticles have charge — e and the lattice 
period of the corresponding QBM model is always one. In the Kane-Fisher model, however, Luttinger-liquid correlation 
g can take any value, and the potential period changes according to g. Another important difference is that in the 
Furusaki-Matveev model, the single particle energy level spacing AE of the dot is much less than T. Therefore, the 
resonance states in the quantum dot form a continuum. On the contrary, AE T in the Kane-Fisher model and 
the double-barrier structure has only one resonance state. Consequently, the on-resonance conductance in the Kane- 
Fisher model is given by the maximum perfect value for a single channel, G* = e 2 /h, whereas in the Furusaki-Matveev 
model, the two QPCs behave like two independent resistors in series, yielding G* — e 2 /2ft,. The difference between 
the two models becomes even clearer when electrons carry spin. In Kane-Fisher model, the spin excitation in the dot 
has a gap, which is the level spacing AE, but the spin is gapless in the Furusaki-Matveev model since AE <C T. 

So far, we have considered only ID QBM models, where there is only one dynamic variable, namely, Q. If one needs 
more than one variable for analysis, the above arguments may be extended to higher dimensional QBM models. Such 
is the case when the spin or transverse degrees of freedom are included into the above two-lead quantum dot model, 
or when the quantum dot is coupled to multiple leads. As far as low-energy physics is concerned, each ID electron 
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channel may be considered equivalent to a lead connected to the dot via a QPC with one transmitting ID channel. 
For example, if there are two spin-channels in the two-lead model, it is equivalent to a four-lead model without spin. 

We may define the number of electrons Qi (i = 1 • • ■ N) for each individual ID channel, where N is the number 
of channels. The set of all variables, (Qi, Q2, • • • Qn), can be thought of as the position vector of an iV-dimensional 
quantum Brownian particle. In the weak coupling limit (w n /r c 3> T, e 2 /C), all Q's are strongly quantized at integers. 
The Brownian particle thus stays most of the time on iV-dimensional cubic lattice sites with unit lattice constant 
(a = 1). As in the above analysis, the number of electrons in the dot, n = const. — JV Qi, may be integrated out at 
low temperatures (T <§; e 2 /C). In this process, the dimension of the resulting QBM model is reduced by one and the 
cubic lattice is projected to the N — 1 dimensions. 

As the two-lead spinless problem of resonant tunneling is mapped to the two-channel Kondo model, so does the 
A-channel resonant tunneling problem to the A^-channcl Kondo model. The mapping will be performed later in two 



steps. Firstly, in Sec. IV, we will map the AT-lead resonant tunneling problem to a QBM model. Then, in Sec. VIII, 
the A^-channel Kondo model will be mapped to the same QBM model. As in the two-lead case, the resonance line 
shape follows a universal scaling function, of which form may be in principle obtained using the properties of the 
iV-channel Kondo fixed point. 



III. RESONANT TUNNELING OF SPINLESS ELECTRONS 



In this section, we will formally construct the QBM action for the two-lead spinless quantum dot model discussed 
in the previous section. Experimentally, one may completely polarize the spin by means of a strong magnetic field. 
In order to avoid quantum Hall effect, we assume that the magnetic field is parallel to the 2DEG plane. As in Ref. 
[| we will assume AE C T< e 2 /C, so quantum coherence is thermally destroyed as electrons travel across the dot. 
This allows us to treat the two QPCs independently. We also suppose that the QjPCs transmit only one ID channel. 
Following the standard bosonization method,c3 we write the Euclidean action asQ 



S = So + S v + Sc 



(3.1) 



where the three terms describe the kinetic energy, backscattering, and the Coulomb energy, respectively. In our units, 

h = Hb = 1. 

The kinetic term is explicitly given by 



(v F d x 8 a ) 



(3.2) 



where 9 a is the boson field of the a-th lead (a = 1,2) and vf is the Fermi velocity of the electrons. In accordance 
with Ref. [| each QPC has its own coordinate system. In particular, we set x — at the center of each QPC and 
assume x < in the dot and x > in the leads, Backscattering is assumed to occur right at x = 0. Since the density 
of electrons is given in terms of boson fields bjO'ESO 



Pa(x) 



d x a (x) 



2n 



(3.3) 



the total number of electrons in the a-th lead is given by 



dx p a {x) 

-0«(O)/27T 



(3.4) 
(3.5) 



up to a constant. 

The backscattering term is given by 



2 

<E 

a=l 



dr r 

T, 



e «.(*=°) +c .c. 



(3.6) 



where v is a dimensionless backscattering amplitude. This describes an electron backscattered from momentum fcj? to 
— kp and vice versa (2fcF-backscattering). In real syste ms, h igher harmonics that describe multiple electron scattering 
(4fci?-backscattering etc.) are also present [See Eq. ( 2.1 C| ) ] . As will be shown in Sec. VI, however, they are less 
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relevant in the renormalization group sense, and it suffices to keep only the single-electron backscattering term in 
order to discuss low-energy physics. We have also assumed that the two QPCs are identical and thus have the same 
backscattering amplitude. The effect of asymmetric QPCs will be discussed later at the end of this section. In general, 
v is a complex number, but without loss of generality it can be always made real by shifting 9 a by a constant. 

Finally, the Coulomb energy term can also be written in terms of 9 a . Since the number of electrons in the dot is 
n = — (Qi + Q2) up to a constant, we may write 



Sc=zc Idr 



9 1 (x = 0)+9 2 (x = 0) 

no 



2tt 



(3.7) 



where no is proportional to the chemical potential of the dot and is controlled by the gate voltage. Note that no is 
not integer in general. 

Now, we may derive the QBM action by integrating out all degrees of freedom but Q\ and Qi. The resulting 
effective action is written as 



S off [Qi,Q 2 ] = \ jdw M [IQiHI 2 + |Q 2 H| 2 ] 
[ dr 

- 2v J — [cos2ttQi(t) + cos27t<3 2 (t)] 



e 2 



+ jcJ dT {_[Ql(r) + Qa(r)] ~ no}2 - (3 - 8) 

We have replaced the Matsubara frequency sum to an integral, which is a valid approximation at low temperatures. 
Identifying (Qi,Q 2 ) with the position vector of a Brownian particle, the above action describes a two-dimensional 
(2D) QBM model. The second term above is now mapped to a periodic potential. Clearly, it has a square lattice 
symmetry with the period a = 1. The minima of the periodic potential are explicitly depicted by circles and squares 
in Q1-Q2 space in Fig. |^. For weak backscattering, one may treat it perturbatively in small v. In the opposite limit 
(large v), the problem is better described in terms of tunneling between potential minima with a small tunneling 
amplitude t. This duality will be discussed in more detail in Sec. fV|. 

If e 2 /C 3> T,v/t c , then n = — (Q1+Q2) can be safely integrated out. In terms of the new variable r = {Q%—Qi)/v^, 
the final effective action is written as 

S eS [r] = l J dLo\Lo\\r{Lo)\ 2 

f dr 

— 4vcos7mo / — cos V^irrM. (3-9) 

J T c 

Note that the above action has the same form as that of the ID QBM model in Sec. [5], except for an additional^ 
dependent factor in front of the periodic potential and a longer period a = y/2. According to Fisher and Zwerger,li3 if 
a > 1, the mobility /i vanishes at T = for any arbitrarily small amplitude of the periodic potential. In other words, 
the Brownian particle gets completely trapped in a potential minimum. In the original quantum dot problem, this 
means that electron tunneling through the dot is completely suppressed at T = 0, unless cos7mo = 0. This is the 
familiar Coulomb blockade. 

If the gate voltage is tuned so that no is exactly a half integer, the periodic potential is identically zero. This is 
the resonance condition where two charge states n = no ± 1/2 become degenerate. In the QBM description, this 
corresponds to a degeneracy between two different sets of lattice sites in the Q1-Q2 space. The degenerate lattice sites 
in Q1-Q2 space forms a "corrugated ID lattice" as in Fig. ||. Note that the period in the r-direction has now become 
a = 1/V2. The motion in the n-direction gets renormalized down to zero at T = 0. Since a < 1, the Brownian particle 
now diffuses freely completely unaffected by the periodic potential at T — 0. This implies that Coulomb blockade is 
lifted and electron current is allowed to flow across the dot. 

We want to make a few comments here. Firstly, the resonance fixed point of the two-lead spinless model is the 
same as the trivial backscattering- free fixed point (v = 0). In general, this is not true in higher dimensional QBM 
models, in which the resonance fixed point is usually inaccessible using a perturbative method. In this case, however, 
the on-resonance conductance is simply given by the perfect conductance calculated in the absence of backscattering. 
From the series conductance of two perfect conductors, we get G — e 2 /2h. 

Secondly, the resonance line shape follows a s cali ng function. If we define <5no as the deviation of no from the 



nearest half integer, one may Taylor expand Eq. (3.E) for small 5no to get 



S v w ±4ttv Sn I — cos V2nr + 0{8nl). (3.10) 

7~r 
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As will be shown in Sec. [vj, the scaling dimension of Sno is given by 1/a 2 . Substituting a = y/2, we find that the 
conductance as a function of Suq and temperature has the scaling form 

G(5n ,T) = G(^). (3.11) 



This is consistent with the result in Rcf. § [See Eq. dglj) ]. Clearly, G(0) = e 2 /2h and G(oo) = 0. As will be shown 



in Sec. VIII , if we map this model to a two-channel Kondo problem, 5no is mapped to the magnetic field TL that is 
locally applied to the magnetic impurity. Since the scaling dimension of Tt is 1/2 [Ref. [33|| , it is consistent with the 
above result. 

Lastly, we briefly discuss the effect of asymmetric QPCs. If v\ ^ V2, there is no longer a symmetry between 
n = tiq ± 1/2 even if hq is a half integer. In fact, the periodic potential never vanishes at any hq. Consequently, 
even the peak value of the conductance resonance vanishes as T — *■ 0. In this sense, there is no true resonance for 
asymmetric QPCs. In the renormalization group point of view, the resonance fixed point is unstable with respect to 
the asymmetry. However, even if the QPCs are not perfectly symmetric, if \v\ — vi\jr c -C T, the flow will not have 
crossed over to the asymmetric fixed point and the system is still governed by the resonance fixed point. 



IV. MULTI-LEAD RESONANT TUNNELING MODELS 



In this section, we will generalize the above analysis to a general TV-lead case and show that it is mapped to a 
(N — l)-dimensional QBM model. 

The first example is the three-lead model shown in Fig. |^(a). As in the two-lead case, we define the total number of 
electrons in each lead Q a (a = 1, 2, 3) and treat them as vector components of the three-dimensional (3D) position of 
a Brownian particle. Charge conservation law again asserts that n = — (Qi + Q2 + Q 3) u p to a constant. Assuming the 



three QPCs are symmetric, the total QBM action looks similar to the one in Eq. (3.8) except for the third variable 
Q3. Again, in the strong backscattering limit (v/t c 3> T, e 2 /G), Q a are strongly quantized in order to minimize 
the potential energy. Since the potential has a cubic lattice symmetry with lattice constant a = 1, the minima also 
form a cubic lattice [Fig. 0(b)]. The Brownian particle propagates by hopping between these lattice sites. At low 
temperatures (T <C e 2 /C), n is integrated out and takes an integer value that minimizes the Coulomb energy. In 
terms of Q a , this condition is 



Q1+Q2 + Q 3 



integer const., 



(4.1) 



which defines a plane perpendicular to the (111) direction. Fig. ||(b) shows two examples of such planes viewed from 
the (111) direction. On any plane, the subset of lattice sites clearly form a triangular lattice with lattice constant 
a = V2. Off resonance, n is allowed to take only one value, and the Brownian particle stays on only one plane. 
However, when uq is a half integer (on resonance), the system is allowed to hop between two planes corresponding 
to n = no ± 1/2. As sketched in Fig. ||(b), they form a "corrugated honeycomb lattice" composed of two parallel 
triangular sublattice planes. 

We may generalize the above analysis to an iV-lead system. In the strong backscattering limit, the coordinates 
(Qi, Q2, ■ ■ ■ , Qn) form an iV-dimensional generalization of a cubic lattice with a = 1. When e 2 /C 3> T, depending 
on whether off or on resonance, the Brownian particle is constrained within one or two (N — l)-dimensional lattices 
on "hyper-planes" perpendicular to the (11- • -1) direction. It is straightforward to show that the lattice on each 
hyper-plane is an (N — l)-dimensional symmorphi(S3 close-packed lattice with a = \J~2~. We will call it an ll (N — 1)- 
dimensional triangular lattice". On resonance (no = half integer), two adjacent hyper-planes are degenerate and the 
two (iV— l)-dimensional triangular sublattices are combined to form an "(N— l)-dimensional honeycomb lattice" . For 
example, if ./V = 4, a 3D triangular lattice (off resonance) is, in fact, a face-centered cubic (FCC) lattice,E3 and a 3D 
honeycomb lattice (on resonance) is a diamond lattice. Note that a diamond lattice is made of two FCC sublattices. 

Now we present a formal derivation of the QBM action from the ./V-lead resonant tunneling model. Generalizing 
Eq. (3.8), we may write 



N M r 

SM{Qa}}=J2{2 J 



dr 

dco \u>\\Q a (uj)\ 2 - v I — cos27rQ a (r) 



2C ldT 



N 



£Qa(T) 



n 



(4.2) 
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For convenience, we perform change of variables 



v 



= 



Qa 



(i=l..-JV-l), 



(4.3) 



where is an orthogonal matrix which is defined in such a way that the last component of the left-hand side of the 
above equation is given by 



N 

After v is integrated out, the effective action becomes 

JV-l 



(4.4) 



S cS [{n}} = i J duj\uj\Y^ \n 



N 



i=l 

(It 



cos 2ir 



7V-1 



no 
N 



(4.5) 



which is an (N — l)-dimensional QBM action. The minima of the periodic potential in {j-j} space form an (N — 1)- 
dimensional triangular lattice unless no is a half integer, in which case they form an (N — l)-dimensional honeycomb 
lattice. The zero-temperature phase diagrams of these lattice models will be derived in Sec. VII. 



V. RESONANT TUNNELING IN A LUTTINGER LIQUID 



A similar analysis as above may be performed for resonant tunneling of a ID Luttinger liquid in a double-barrier 
structure. Let us first discuss the spinless case. Following Ref. [ll| we suppose that there are (5-function barriers at 
x = and x — d. The Euclidean action in the bosonized form is given by 



S 



8irgv s 



drdx [(d T 9) 2 + (v s d x 9) 2 ] + eV G / dr [6(d) - 0(0)] + v — {cos[0(O)] + cos[0(d)]} 



d,T 



(5.1) 



where v s is the sound velocity, g is the Luttinger-liquid correlation, and Vg is the gate voltage near the isolated region 
< x < d. If the electron-electron interaction is repulsive (attractive), the Luttinger-liquid correlation parameter 
satisfies g < 1 (g > 1). After 9(x) is integrated out except at x = and d, we get the effective action 



S e s[Q,n] = - / dui \w\ 

gJ 



1 



|n(w)| 2 



'\d/v s 



4 \ — g— \u\d/v s 



bVq J dr n(r) + v J — cos2ttQ(t) cos7rn(r), (5.2) 



where 



Q = 



6(0) + fl(rf) 
4tt 

9(d) - 6(0) 
2tt 



(5.3) 
(5.4) 



Q is the number of electrons that has transferred from one side of the double-barrier to the other and n is the number 
of electrons between the barriers. 

Since there is an energy gap AE = v s /gd for n(u>) as ui — > 0, we may safely integrate it out for |cj| <C AE. The low 
frequency part of the action is then written as 



Scs[r] — — J cLa \uj\\r(uj)\ 2 + v cos irriQ J — cos 2tt y/gr(r), 



(5.5) 



where r = Q/y/g and no = eVGgd/v s . 

The above action closely resembles that of the two-lead quantum dot model in Eq. (3.9). However, there are a 
few important differences. First of all, the potential period is 1/-Jg in the Luttinger liquid, which varies with the 
electron-electron interaction, whereas it is fixed at y2 in the quantum dot model. Another important difference is 
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that the level spacing AE is finite in the Luttinger liquid model. Thus, at low temperatures (T <C AE), the thermal 
coherence length Lt = Tiv s /kBT is much greater than the barrier separation d. In other words, when electrons tunnel 
through the double-barrier, they fully retain their phase information. Therefore, a single boson field 9{x) has been 
used for both barriers, whereas in the quantum dot model, each QPC was described by an independent field 9 a {x). 
Consequently, the two models are not the same even in the non- interacting limit of the Luttinger liquid (g = 1). In 
that limit, the period is one for the Luttinger liquid model and the on-resonance conductance is given by G* = e 2 /h, 
which differs from that of the two-lead quantum dot model by a factor of 2. The resonance line shape also has a 
different scaling form, 



G(Sn Q ,T) 



\ T i- 



(5.6) 



Interestingly, the Luttinger liquid model yields the same on-resonance conductance and resonance scaling function as 
the two-lead quantum dot model, if g = 1/2. 

A similar analysis is possible when there is spin degree of freedom. Now we need two boson fields 0j and &i for up 
and down spin states. It is convenient to define charge and spin fields 



h-h- 



The effective action is given by 
S e ff[{Q a },{ n <x}]= — / duj 



\Qa(u) 



|n Q (w)| 5 



1 



-\uj\d/v s 



4 1 



-\uj\d/v s 



(It 



(5.7) 
(5.8) 



el/ G n,(T)+H^(r)]+S„, (5.9) 



where TL is the magnetic field applied only between the barriers and S v contains all allowed backscattering processes. g p 
and g (j are Luttinger-liquid correlations for the charge and spin sectors, respectively. If there is SU(2) spin symmetry, 
g a = 2. In the non- interacting limit, g p — g a = 2. Physically, Q a means the total spin that has transferred across the 
barriers and n a means the total spin between the barriers (both in units of h/2). Q p and n p have similar meanings 
except that they are defined for number of electrons. Note that both n p and n a have an energy gap, whereas in the 
quantum dot model, only n p has a gap. 

After both n p and n a are integrated out, the above effective action describes a 2D QBM in terms of Q p and Q a . 
Particularly in the vicinity of g p = 1 and g a — 3, Kane and-Edsher have found a fixed point similar to the intermediate 
fixed point of the QBM model on a 2D honeycomb latticeJl3 There, the most relevant backscattering terms are 



r a] T 

S v = I — [v e cos ttQ p cos ttQo- + v\ cos 2irQ p + V2 sin 2irQ p ] , 

J Tc 



(5.10) 



where all u's are real-valued. Physically, the first term describes backscattering of a single electron and the other two 
backscattering of a spin-singlet pair. Defining rescaled variables 



the low-frequency effective action is written as 



Qa/y/g^, 



(5.11) 



1 f f di~ 

S c s[r P , r a ] — - j du)\u>\ [|r p (w)| 2 + |r CT (w)| 2 ] + / — (y e cosTTy^r p cosiry/g^ra + v\ cos 2-k ^fg~ p r p + v 2 sm2Tr y /g p ~r p ) . 

(5.12) 



if 9<?/g P 



3, V e 



2vi, and v-2. = 0, the above action has a honeycomb lattice symmetry in the r p -r a plane. For 



9<r/9p 7^ 3, the lattice is stretched or compressed. The fixed point of this "distorted honeycomb lattice" will be 
discussed in Sec. VII E . 



VI. GENERAL FORMALISM OF THE QBM MODEL 

In this section, we will present a general formalism of the QBM model with a periodic potential and its renormal- 
ization group analysis. Specific lattice models will be discussed in the following section. In any arbitrary dimensions 
a generic QBM action is written as 
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S — So + S v , 



(6.1) 



where 

So = \ J du MeM^rHI 2 , (6.2) 

& = - / ^^voe^-'W. (6.3) 

c G 

The above two terms represent the dissipative kinetic energy and the periodic potential, respectively. r(r) is the 

trajectory of the Brownian particle as a function of imaginary time r and vg is a dimensionless amplitude of the 

potential at a reciprocal lattice vector G. In Eq. ( |6.2[ ), high frequency fluctuations are exponentially suppressed by 
t c . The potential minima form a lattice {R}, where 

G R = integer, (6.4) 



for any G and R. The dissipation is proportional to the coefficient in the righthand side of Eq. (6.2), but it may be 
absorbed into G by rescaling r and G appropriately. As dissipation gets stronger, |R| increases and |G| decreases. 

Now let us discuss the linear response of the Brown particle to a uniform external force F. Firstly, we add the 
following term to the total action: 

S F = -J drF-r(r). (6.5) 

Then the mobility ynj is defined as the ratio of the average velocity of the Brownian particle to the external force, in 
the linear response regime. Following standard linear response theory, it is written as 

= lim — !- / du' ajui'(ri(u)rj(— ui 1 )), (6.7) 
uj^o M J 

(6.8) 

In the last line we assumed that the mobility is isotropic. Note that it is normalized in such a way that fi — 1 ("perfect 
mobility"), if Vg =0. 

In the weak potential limit (vg *C 1), it is natural to construct a perturbation theory in small vg- In order to 
study low-temperature properties, we perform a renormalization group analysis analogous to the ones in Ref. EM The 
general scheme of the renormalization group analysis is described in Appendix A. Up to the leading order in vg, the 
renormalization group flow equation is given by 

= (1 - |G|> G . (6.9) 

Clearly, all vg are irrelevant if the shortest reciprocal lattice vector satisfies |Go| 2 < 1. 

On the other hand, if vg S> 1, the Brownian particle is more likely to be localized at a minimum R of the potential. 
In that case, transport is dominated by tunneling between potential minima, and the partition function may be 
expanded in powers of the fugacity of tunneling events. In order to describe tunneling, it is more convenient to write 
the action in terms of k(r), which is the trajectory of the particle in the momentum space. The new action, which is 



dual to Eq. (6.1), is given by 



where 



S = S + S t , (6.10) 

5 = i [du>\w\eW°\k(u>)\ 2 , (6.11) 
S t = - [ -Vi R e«- k W. (6.12) 
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The two terms represent the kinetic energy and the tunneling, respectively. £r is the tunneling amplitude between 
two lattice sites connected by a displacement vector R. The mobility is written as 

fiij = Sij — lim - — - ( dio uiu! 1 ' (ki{uj)kj{—uj')) . (6.13) 
u^o \uj\ J 

Since the dual action has exactly same form as the original action, it is straightforward to obtain the renormalization 
group flow equation 

^: = (l-|R| 2 )i R . (6.14) 

From the above flow equations one can easily find two trivial fixed points: the localized (£r, = 0) and the free 
(i>g = 0) fixed points. The mobility at the fixed points is simply given by /i = and fj, = 1, respectively. Other fixed 
points will be obtained and analyzed in the next section. 



VII. LATTICE MODELS AND PHASE DIAGRAMS 



In this section, we study the properties of the QBM model with various lattice symmetries. Using a renormalization 
group analysis, we will derive the phase diagram of each lattice. 

A lattice is either symmorphic or non-symmorphic depending on the crystallographic space group of the coordinate 
transformations that leave the system invariant. If the space group is a semi-direct product of— a pure translation 
group and a pure rotation group, it is called symmorphic. Otherwise, it is non-symmorphic.Ej For example, all 
Bravais lattices are symmorphic, and the honeycomb, Kagome, and diamond lattices are non-symmorphic. One of 
the major results of this paper is the proof that there is a stable intermediate fixed point only for non-symmorphic 
lattices. The first two subsections below deal with symmorphic lattices and the following three with non-symmorphic 
lattices. 

A. One-dimensional lattice 

The simplest example of a symmorphic lattice is a ID lattice. According to Eq. (^J) the lattice sites in the direct 
and the reciprocal lattice are given by 

R = mR , (7.1) 
G = n/Ro, (7.2) 



where Rq is the lattice constant, and m and n are integers. It follows from Eqs. (|6.9| ) and ( |6.14| ) that 



^ = (l-n 2 ARo> G , (n^O), (7.3) 

^ = (1 - m 2 Rl)t R , (m^O). (7.4) 

From the above equations, it is obvious that if i?o < 1, all vq are irrelevant, whereas is relevant at least for 
rn = ±1. Now suppose that the bare value of tn(m — ±1) is perturbatively small. If the system is scaled down to 
lower energies under successive renormalization group transformations, tn will keep growing until eventually it departs 
from the perturbative regime. Presumably, the system will continuously flow towards the opposite limit where vq 
is small. Under further renormalization group transformations, vg will decrease until it vanish at T = 0. In more 
formal words, the localized fixed point (tn — 0) is unstable and the free fixed point (vq = 0) is stable. At T = 0, the 
Brownian particle diffuses completely unaffected by the periodic potential. 

In contrast, if Rq > 1, all ta are irrelevant whereas vg is relevant at least for n — ±1. Now, exactly the opposite 
happens under renormalization group transformations. In short, the free fixed point {vq = 0) is unstable and the 
localized fixed point (tji = 0) is stable. The zero-temperature phase is now localized and the particle gets completely 
trapped in a minimum of the potential at T = 0. 

The physical interpretation of these is as follows. For R Q < 1, the lattice sites are close enough for a given dissipation 
strength — or the dissipation is weak enough for a given lattice constant — that the Brownian particle tunnels easily 
at low temperatures. If Rq > 1, however, lattice sites are far apart or the dissipation is too strong that the particle 
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can hardly tunnel. In fact, the above result is a reproduction of the result of a ID QBM model in Ref. [L6| and 
of a single-barrier spinless Luttinger liquid model in Ref. [l8|. In the latter model, the Luttinger-liquid correlation 
parameter g is related to the lattice constant by 

g = 1/Rl (7.5) 

The above free and localized phases are equivalent to the fully coherently transmitting phase of the Luttinger liquid 
model for attractive interactions (g > 1, Rq < 1) and the insulating phase for repulsive interactions (g < 1, Rq > 1), 
respectively. 

The above result may be summarized in terms of the zero-temperature universal mobility 

r o if Go < i, , , 

M \ 1 if G > 1. { '- 0> 

Note that the zero-temperature phase is completely independent of v or t. The phase transition is solely controlled 
by R . 

In order to determine the zero-temperature phase, it is sufficient to consider only the shortest lattice vectors Ro 
and Go, because they are the most relevant operators in the respective perturbative regime. Throughout this paper, 
we will ignore all the other less relevant lattice vectors. However, one can easily include less relevant operators when 
it is necessary. 

Without many modifications, the above analysis may be applied to square lattices, cubic lattices, and their higher- 
dimensional generalization. In the vicinity of the phase boundary |Go| 2 = 1, each vector component of r or k 
gets completely decoupled from the others at T = 0, since all the coupling terms such as e i27r( - ri+r ^ are irrelevant. 
Therefore, the square lattice and its higher dimensional generalizations can be treated as a set of independent ID 
lattices. 



B. Triangular lattice and its D-dimensional generalization 



We argued in Sec. IV that the off-resonance multi-lead quantum dot model is mapped to the QBM model on a 
multi-dimensional triangular lattice. In this subsection, we will first analyze the 2D triangular lattice model, and then 
generalize it to higher dimensions. 

The reciprocal lattice of a triangular lattice is also a triangular lattice. We will assume that there is a six-fold 
symmetry and the potential amplitudes are positive for the six shortest reciprocal lattice vectors Go, i.e., 

lid, = V — const. > 0. (7.7) 

From the six-fold symmetry, it also follows that 

^R = t = const. (7.8) 

It is straightforward to see that 

|G | 2 |Ro| 2 = | (7.9) 

A simple renormalization group analysis as in the previous subsection shows that if | Go | 2 > 4/3 and |Ro| 2 < 1, the 
free fixed point is stable and the localized fixed point in unstable. If |Go| 2 < 1 and |Ro| 2 > 4/3, the localized fixed 
point is stable and the free fixed point is unstable. However, if 1 < |Go| 2 , |Ro| 2 < 4/3, both fixed points are stable 
and there is at least one unstable fixed point between them. Physically, this means that there are at least two distinct 
phases at T — 0: free and localized. Although it is very unlikely that there are more fixed points, we currently do 
not have an evidence to completely rule out the possibility. Assuming there are only two stable fixed points, the 
zero-temperature phase depends on whether the bare value of v is above or below the unstable intermediate fixed 
point. As v moves from one side of the unstable fixed point to the other, a first order phase transition occurs at T = 
and the mobility jumps between zero and one. 

For small v, the intermediate fixed point may be perturbatively accessed near | Go | 2 = 1- More specifically, we 
introduce a small parameter e > and perform a perturbative renormalization group analysis at | Go | 2 = 1 + e. Since 
details of the analysis will be given in Appendix A, here we will just state the result. The renormalization group flow 
equation is given by 
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dv 9 . 

— = -ev + 2v 2 . (7.10) 



Solving dv/d£ = 0, we get the intermediate fixed point 



|. (7.H) 



The scaling dimension of the most relevant operator at the fixed point is 

A = l-c. (7.12) 
Using Eq. we may also compute the universal fixed point mobility 

M* = 1 " ^e 2 . (7.13) 

A similar perturbative analysis was used in Ref. [l8]. 

The analysis for small t near |Ro| 2 = 1 is very similar. In fact, since both the direct and reciprocal lattices are 
triangular, one may directly use the duality. Namely, the small v limit is mapped exactly to the small t limit if the 
following substitutions are made: 

r^k, (7.14) 

G -> R, (7.15) 

u->t. (7.16) 



The mobi lity /j, is then mapped to 1 — \i. At |R | 2 = 1 + e, the results can be directly read off from Eq. ( 7.10 ) through 
Eq. (|7.13[). Explicitly, they are 



a— ^ 

t* = |, (7.18) 

A=l-e, (7.19) 
3tt 2 

M* = ^e 2 . (7.20) 

Note that there is a self-dual point at |Go| 2 = |Ro| 2 = v3/2 where the dual models become identical. The mobility 
at the self-dual point is obviously given by 

= 1 " A4i = 1/2. (7.21) 

It is very hard to calculate the intermediate fixed point away from the above two perturbative limits and the self- 
dual point. However, assuming that the fixed-point lines on both limits smoothly join together at the self-dual point, 
one may complete the phase diagram as shown in Fig. |^(a). The top and bottom lines of the figure depicts the free 
(v = 0) and the localized (t — 0) fixed points, respectively. The intermediate fixed-point line is also drawn along with 
the renormalization group flows. The above results may be summarized in terms of the fixed-point mobility 

I if|G | 2 <l, 

fi* = { or 1 if 1 < |G | 2 < 4/3, (7.22) 
if |G | 2 >4/3. 

Unlike in the ID lattice problem, v and t control a first order phase transition, if 1 < | Go | 2 < 4/3. 

Now, let us generalize the above analysis to the D-dimensional triangular lattice. As discussed in Sec. IV, a D- 
dimensional triangular lattice is defined as a subset of the lattice sites of a (D + l)-dimensional cubic lattice, that are 
constrained to a hyper-plane perpendicular to the (11 • • • 1) direction. Its lattice sites are equivalent to the minima 
of the off-resonance backscattering action of the (D + l)-lead quantum dot model [Eq. Q4.5Q with n ^ half integer], 
except for that |Ro| is no more fixed at Specifically, it is a conventional triangular lattice for D — 2 and an 
FCC for D — 3. It is hard to visualize an object in dimensions greater than three, but it is helpful to illustrate 
the D-dimensional triangular lattice using an induction argument. For example, we can construct an FCC from a 
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2D triangular lattice in the following way: A triangular lattice is made of periodically placed equilateral triangles. 
Likewise, an FCC is made of equilateral tetrahedrons. A tetrahedron in an FCC consists of a triangle at the bottom 
and the fourth vertex placed above the triangle plane. Due to periodicity, this last vertex also belongs to a triangle 
in another plane parallel to the first. In other words, an FCC is obtained by stacking up planes of triangular lattices. 
However, it has to be done in a certain way that the resulting lattice is a symmorphic close-packed structure.^ In 
general, a ZJ-dimensional triangular lattice is obtained by stacking up (D — l)-dimensional triangular lattice hyper- 
planes. 

Now we may perform a similar renormalization group analysis as above. If D > 3, however, the reciprocal lattice 
{G} is no longer the same as the direct lattice {R} and we cannot use duality. (For example, the reciprocal of an 
FCC is a BCC.) From a geometrical consideration, we find 

2 D 

|G | 2 |Ro| 2 = (7.23) 

If 1 < | Go | 2 < 2D/(D + 1), there is an unstable fixed point that separates the localized and the free phases. The 
phase diagram is schematically drawn in Fig. ^|(b). We may also perform an e-expansion analysis analogous to the 
above. The results are listed in Table [j] at | Go | 2 = 1 + e in the small v limit. In the opposite limit of small t, the 
results are given in Table H at |R | 2 = 1 + e. 

Before we move on to non-symmorphic lattices, we would like to discuss briefly that the symmetry of the lattice 
is crucial on determining the phase diagram. For any arbitrary Bravais lattice, |Gq|| Rq| > 1, and obviously, at least 
one of | Go | and |Rq| is greater than one. Therefore, it follows from Eqs. (3.£) and ( |6.14 ) that at least one trivial 



fixed point is stable for any arbitrary lattice constant. On the contrary, for a non-Bravais lattice, | Go | and |Rq| 
may be less than one at the same time, in which case both trivial fixed points are unstable. Then, there must be 
a stable intermediate fixed point between them. It is indeed the case for the non-symmorphic lattices such as the 
honeycomb and the diamond lattice. A common feature of all non-symmorphic lattices is that they are made of Bravais 
sublattices and there is a geometrical symmetry between the sublattices. Since the asymmetry between sublattices 
is always relevant, even if the symmetry is slightly broken, one sublattice will dominate the others after repetitive 
renormalization group transformations. Then, the problem reduces to that of the dominating sublattice, which is a 
symmorphic lattice. Therefore, it is the extra symmetry between sublattices that prevents this "renormalization of 
basis", and non-symmorphicity is a necessary condition for the existence of a stable intermediate fixed point. 



C. Honeycomb lattice and its D-dimensional generalization 



In this and the following two subsections, we will study non-symmorphic lattice models. The first example that will 
be discussed in this subsection is the Z?-dimensional honeycomb lattice. Following the framework of the the previous 
subsection, we will first consider the conventional 2D honeycomb lattice, and then generalize it to higher dimensions. 

A honeycomb lattice is made of two triangular sublattices. As depicted in Fig. g, the six shortest lattice vectors 
Ro € {±R a } (a = 1,2,3) always connect one sublattice to the other. Note that Ro € {R a } on A sublattice and 
Ro £ {— R a } on B sublattice. The reciprocal lattice {G} is a triangular lattice and the six shortest reciprocal lattice 
vectors are drawn in Fig. ^(a) as Go 6 {±G a } (a = 1, 2, 3). It is straightforward to see that 

|Go| 2 |Ro| 2 ^. (7.24) 



From Eqs. (6J)) and ( 6.14 ), it follows that for 4/9 < | Go | 2 < 1, both the free and localized fixed points are unstable. 
It implies that there must be at least one stable fixed point between the two. Again, we will assume that there is 
only one stable fixed point. At the new stable fixed point, the mobility takes a universal value between zero and one 
(0 < < 1). 

Now let us perform an e-expansion analysis similar to the one in the previous subsection. Computing the structure 
factor, the potential amplitudes are given by 

VGo = „(1 + e ^G .R 0) = ue ±<ir/3 i (? 25) 

where we have chosen the origin of the 2D space to lie on an A-sublattice site. The sign in the exponent is (+) if 
Go € {G a } and it is (— ) if Go € {— G a }. Using the above equation, the periodic potential part of the action may be 
written as 
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dr 



i27rGo -r 



Gn 



2v J ^£cos(27rG a .r-| 

a— 1 



(7.26) 
(7.27) 



Since a trip around a triangular plaquette in the reciprocal lattice space gives rise to a phase ±7r [Fig. ^(b)], this 
problem may be also viewed as a tight-binding representation of a triangular lattice in which a flux it is threaded in 
each plaquette with alternating sign. 

The e-expansion analysis in the small v lim it near | Go | 2 = 1 is essentially the same as that of a triangular lattice 
except for the extra factor e ±l7r / 3 in Eq. ( 7.25| ). The results are summarized in Table |. In the opposite limit of small 
t near |Ro| 2 = 1, however, the situation is more complicated due to the two sublattices. When many tunneling events 
occur sequentially, Ro must be alternately chosen from {R a } and {— R a }. Technically, this is automatically taken 
care of by using the Pauli matrices ± ia y . The tunneling part of the action may be written as 



S t = -t 



i27rR„-k 



(7.28) 



The alternation between {R a } and {— R a } is ensured by the identity (cr*) 2 
is given in Appendix A. The results may be summarized as follows: 

dt „ a 

di =et - M > 



t* 



A = 1 + 2e, 

* 2 

/i = 7T e. 



0. Details of the e-expansion calculation 

(7.29) 

(7.30) 

(7.31) 
(7.32) 



Fig. |](a) shows the phase diagram of the honeycomb lattice. Since the stable intermediate fixed-point line smoothly 
connects the localized and the free phases, there is no first order phase transition as in symmorphic lattices. The fixed 
point mobility //* is also a continuous function of | Go | 2 - In general, the calculation of /i* is not easy for an arbitrary 
|G | 2 , but as will be shown in the Sec. VIII, an exact solution may be obtained for the special case where |G | 2 = 2/3, 



using the mapping of the honeycomb lattice model to the Toulouse limit of the three-channel Kondo model. 

Now let us generalize the above analysis to D > 3. First of all, we define a D-dimcnsional honeycomb lattice as a 
subset of the lattice sites of a (D + l)-dimensional cubic lattice, that lie on two adjoining hyper-planes perpendicular 
to the (11- ■ -1) direction. Its lattice sites coincide with the minima of the on-resonance backscattering action of 
the (D + l)-lead quantum dot model [Eq. ( p~5| ) with no = half integer], except for that |Ro| is no longer fixed. For 
example, a "3D honeycomb lattice" is a diamond lattice. As a honeycomb lattice consists of two triangular sublattices 
and a diamond lattice consists of two interpenetrating FCC sublattices, a generic D-dimensional honeycomb lattice 
consists of two interpenetrating Z?-dimensional triangular sublattices. Among the two basis sites of a D-dimcnsional 
honeycomb lattice, if one is placed at a vertex of a D-dimensional equilateral (D + l)-hedron, the other is at the 
center. 



Generalizing Eqs. (7.24) through (7.28), we get 



|Ro| 2 |G | 2 = 



D 



D + l 



D + l 



/dr ( 7r \ 

-I Yi, cos (Z^Ga ■ r - _j 



St = -t (^-e^ R ^ + H.c. 

J Tc a=l \ 



(7.33) 
(7.34) 
(7.35) 

(7.36) 



In terms of the matrix defined in Eq. (4.S), the lattice vectors are explicitly given by 
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G a = c(0- 1 \0- 2 \---,0-l > ), (7.37) 
Ra = -i(0- 1 1 ,0- 2 1 ,--- ) 0^) 

Q 

= £ (c = const.). (7.38) 



Note that Eq. (|7.35|) is identical to the second term in Eq. fl4.5| ) if N = D + 1, no = —1/2, and c = 1. There is a 



stable intermediate fixed point if [-D/(I? + l)] 2 < |Go| 2 < 1- The results of the e-expansion analysis are summarized 
in Tables | and || 

Quite interestingly, the phase diagrams for D < 3 [Fig. ||(a),(b)] and for D > 4 [Fig. ||(c)] have quite different 
structures. Unlike for D = 2 or 3, the fixed-point line for D > 4 approaches the top line from right. Consequently, 
arbitrarily close to the top line, the intermediate fixed points are unstable. As schematically shown in Fig. |8|(c), we 
conjecture that the fixed-point line eventually curves back at a certain critical point £ and gets continuously connected 
to the bottom line. Unfortunately, this critical point is not perturbatively accessible, so a definitive evidence to prove 
(or disprove) the above conjecture is currently not available. In the next subsection, however, we present an argument 
that supports the conjecture, by using an analysis on a distorted diamond lattice. For this exotic model, the direction 
from which the intermediate fixed-point line approaches the top line changes continuously from right to left as the 
lattice is deformed. Thus, if the fixed-point line is to be smoothly connected to the bottom line for D < 3, it is very 
likely that it will also be the case for D > 4 despite the difference, because the different feature (the approaching 
angle) is continuously deformable. In short, the above conjecture is no less reasonable for D > 4 than for D < 3. As 
shown in Fig. ||(c), if 1 < | Go | 2 < £, the free fixed point and one of the two intermediate fixed points are both stable. 
A first order phase transition between the free and the intermediate phases occurs at T = as v or t is tuned. 



D. Distorted diamond lattice 



In this section, we will deform a diamond lattice in a continuous manner to see how the angle with which the 
intermediate fixed line approaches the free fixed-point line changes. A distorted diamond lattice is defined by the 
minima of the potential 



4 

V = -2— ^cos(27rG a -r- |) , (7.39) 

a— 1 



where 



Gi = c(a, 1, 1), 
G 2 = c(a, -1, -1), 
G 3 = c(— a, 1, -1), 

G4 = c(— a, — 1, 1), (c, a = const.). (7.40) 

This lattice is obtained by stretching (or compressing) a diamond lattice in the (100) direction. The regular diamond 
lattice is restored when a = 1. 

One of the shortest displacement vectors connecting the two sublattices is 

Rl =-S& 1 ' 1 )- (7 ' 41) 
Thus, the relation between the lengths of the shortest direct and reciprocal lattice vectors may be written 

ICI-IB.,. . <°' + f< 2 f + 1 > . (7.42, 

If (11 — Vl05)/4 < a 2 < (ll + Vl05)/4, then |G | 2 |Ro| 2 < 1 and both trivial fixed points are unstable. Using a similar 
e-expansion calculation as in the previous subsection, it is straightforward to see that there is a stable intermediate 
fixed point at |Ro| 2 = 1 — e near the localized fixed point (e > 0). The derivation of a flow equation at |Go| 2 = 1 — e 
for small v is more complicated, but it may be eventually written as 

^ = ev-f(a)v 3 . (7.43) 
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The function f(a) is computed numerically. If a < 1.9073, then f(a) > and the solution to the above equation 
exists only if e > 0. The intermediate fixed points are stable and approach the top line from left, as is the case for 
a regular diamond lattice [Fig. ||(b)]. On the other hand, if a > 1.9073, then f(a) < and the solution exists only 
if e > 0. The fixed points are now unstable and approach the top line from right. Consequently, the phase diagram 
resembles that of a four or higher dimensional honeycomb lattice [Fig. ^|(c)]. Since /(a) is a continuous function, the 
angle with which the fixed-point line approaches the top line is also continuous. Therefore, we come to the following 
conclusion: If our conjecture in the previous subsection, that the two trivial fixed-point lines are connected by a 
single intermediate fixed-point line, is correct for D < 3, it is most likely correct for D > 4, too, even though the 
intermediate fixed-point line has a different shape in the two cases. 



E. Distorted honeycomb lattice: resonant tunneling of a spin-1/2 Luttinger liquid 

As argued in Sec. [y], the double-barrier resonant tunneling model of spin-1/2 Luttinger liquid is mapped to the 
distorted honeycomb lattice. Following Ref. [ll| we will briefly discuss the intermediate fixed point of this exotic 
lattice problem. 

Since all v's in Eq. ( |5.12| ) are marginal if g p — 1 and g„ = 3, one can perform a perturbative renormalization group 
analysis at g p = 1 — e p and g a = 3 — e CT , for small positive e's. According to Ref. |l8|, there is a stable intermediate fixed 
point at v* = \/e p (e p + e CT ), v* — (e p + £cr)/4, and v% = 0. The renormalization group flow diagram is schematically 
drawn in Fig. M in the v e -v\ plane, assuming V2 — 0. This intermediate fixed point is a continuation of that of the 
regular honeycomb lattice. However, unlike the regular honeycomb lattice, which is described by only one parameter 
v, there are now three parameters. Since the intermediate fixed point is stable only in one direction, in order to obtain 
resonance, one has to tune not only one, but two parameters, for example, V2 and v e /v\. 



VIII. MULTICHANNEL KONDO PROBLEM 



In this section, we will explicitly derive the mapping between the multichannel Kondo model and the QBM model. 
Then using the mapping, we will exactly calculate the mobility fj,* at the intermediate fixed point-.—. 
In the absence of a magnetic field, the Hamiltonian of the A^-channel Kondo model is given byOEj 

H = H + Hj, (8.1) 

where 

H = ivf^2 / dxq PL9x->Pa S , (8.2) 

a.s 

Hj = 2irv F J2{jzSL p s z a (0) 



■5* [' 



5+ s-(0)+H.c. . (8.3) 



imp 

The first term Hq is the kinetic energy of electrons and Hj is the exchange coupling between the impurity spin <Si mp 
and the sum of electron spins ^ a s a = J2 a i>l s (&ss' /2)^>as' at x = 0. a is the Pauli matrix and a and s are channeLand 
spin indices, respectively. For generality, we assume J z and J± are independent. Following Emery and Kivelson,E3 we 
bosonize the electron operators 

<Pas= J- e~'»", (8.4) 

y'lTTVpTc 

where Vf is the Fermi velocity, r c is the short-time cutoff, and <p as is a boson field that satisfies the following 
commutation relation 

[<p as (x), <j>a's'{x')] = -iir5 aa ,6 ss > sgn(x - x'). (8.5) 
If we separate spin <f>° = tft a ^ — <p a i from charge § p a — 4> a ^ + (f> a ^, we may write s a as 
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si = —d x 



47T 



1 



2ttvft c 

In the boson representation, the Kondo Hamiltonian is written as 

N 



dx 



cr\2 



1 



JV 



o=l 



imp 



+H.c. 



(8.6) 
(8.7) 

(8.8) 
(8.9) 



Since 4>a is free and non- inte racting, we will drop it from now on. 

The first term in Eq. ( |8.9| ) may be eliminated by introducing a unitary transformation 

which rotates the axes of the spin space at x = 0. The Hamiltonian transforms as 

h' = Utfir 1 



E 



8tt 



er\2 



dx(dx<t>l) 



V F 



2r r 



f 1 p 

'-'imp 



H.c. 



(8.10) 

(8.11) 
(8.12) 



Now, we set A = J z /2 so that the second term above vanishes. A similar technique has been used in Ref. |3^ for the 
two-channel Kondo problem and in Ref. [36| for the four-channel problem. 

For later convenience, we make change of variables using the the same orthogonal transformation as in Eq. fl4.3p: 











= 





(i = l---JV-l). 



(8.13) 



The last component of the left-hand side, s = Ea^a/ViV is proportional to the sum of the spins in all channels, 
whereas the other components c6^ (i = 1, • • • , JV — 1) are "spin flavors". In terms of these new variables, the above 
Hamiltonian finally becomes 



H' 



VF 

8tt 



dx 



s\2 



N-l 

i=i 



sf\2 
x^i ) 



Jl 

2t c 



E {5+ pe ^[^( 1 ~ fJ ^ 0S(O)+ ^""°" 0f(O) ] +H, 



a=l 



.14) 



Note that there is a reflection symmetry about J z = 2/JV. More specifically, H' is invariant under (ff — > — </> s and 
J z — * — J z + 4/JV. At the symmetry point J z = 2/JV, the total spin </> s gets completely decoupled andJiecomes 
non-interacting. This is the multichannel generalization of the Toulouse limit of the single-channel modelEj or the 
Emery-Kivelson line of the two-channel modelH Although this special limit is solved exactly for JV < 2, it is not the 
case for JV > 3. As will be shown later, however, there is a stable fixed point in the Toulouse limit, so some of the low- 
energy properties may be understood via a renormalization group analysis. This fixed point is identical to the stable 
intermediate fixed point of the QBM model on an (JV — l)-dimensional honeycomb lattice with | Go | 2 = (JV — 1)/JV. 

Now we may obtain the QBM action using the same technique that we have used for the multi-lead quantum dot 
model. By setting J z = 2/JV and integrating out all degrees of freedom away from x — 0, we get 



SKondo I^W + T / f E 



imp 



(8.15) 



The a bove action has identically the same form as the multi-dimensional honeycomb lattice model in Eqs. (3.11), 
(7.36), and (7.38), with c = 1. The mapping becomes complete if we make the following substitutions: 



N^D + l, 
<f>f(x = 0) 271-fci, 
Ji 2t, 

O-- 1 «-> Rn 



to- 



(8.16) 
(8.17) 
(8.18) 
(8.19) 

(8.20) 
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This mapping may also be understood using the "spin configuration lattice" , just as we have used lattices in 
Qi — Q2 — • ■ ■ — Qn space for the quantum dot problems. The spin configuration lattice lives in a space where each 
axis represents the total spin in channel a, 

Sl=jdxsl{x). (8.21) 

For small coupling J±, the change in will occur in a step-like manner with unit step height (AS*^ = 1). For the 
TV-channel model, therefore the vectors (Sf , 5|, ■ ■ ■) form an iV-dimensional cubic lattice with lattice constant a = 1. 
In this lattice, a spin transfer between channels is represented by tunneling of a fictitious particle between lattice sites. 
The conservation condition of the total electron spin plus the impurity spin gives a constraint in the lattice space. 
Explicitly, the constraint may be written as 

5>*= const. -Sf mp . (8.22) 

a 

This is an equation that defines an (N — l)-dimensional hyper-plane perpendicular to (11 • • • 1) direction. Since 
S? mp = ±1/2, the fictitious particle is allowed to reside on two adjoining hyper-planes that satisfy the above condition. 
For the three-channel Kondo model, the two (111) planes are analogous to those shown in Fig. |](b). In the Toulouse 
limit (J 2 = 2/N), the motion perpendicular to the planes gets renormalized at low temperatures and the two planes 
eventually collapse onto each other at T — 0. As a result, we end up with a new tunneling problem in a single plane. 
In the three-channel case, the resulting lattice is a 2D honeycomb lattice [Fig. f|(b)]. In general, the TV-channel Kondo 
problem in the Toulouse limit is mapped to an (N — l)-dimensional honeycomb lattice-|QBM model. 

A renormalization group analysis similar to that of Anderson, Yuval, and HamanntZI may be performed to obtain 
the flow diagram of the multichannel Kondo problem. Details of the calculations are given in Appendix B and we 
will only present the result here: 

f = 4d-^), (8.23) 

dj± N . N o , 

-j± = J t J ± (l--J z )--Jl (8.24) 

The above equations are exact in J z and perturbative in Jj_. It is also exact in TV, except for the last term in Eq. 
( |3.24 ), in which we kept only the leading order term in 1/N . The renormalization flow diagram is schematically drawn 



in Fig. [T(]. Obviously, the multichannel Kondo fixed point lies on the Toulouse line J z — 2/N. Although it is a highly 
non-trivial strong coupling fixed point, we may use known results of the conformal field theory to exactly compute 
some physical quantities. On the other hand, we may use the mapping to the QBM model and perturbatively access 
the fixed point in the large D limit of the D-dimensional honeycomb lattice model. Therefore, it is possible to directly 
compare the two results, in order to test the mapping and the perturbative analysis of the QBM model in a non-trivial 
way. Specifically, we will calculate below the scaling dimension of the leading irrelevant operator A and the fixed 
point mobility fi* using both methods. 

From a conformal field theory calculation, the scaling dimension of the leading irrelevant operator at the Kondo 
fixed point is given byE3 

(8.25) 



N + 2 

We may also compute the same quantity using an e-expansion technique described in Appendix C. By substituting 



N = D + 1 in Eq. (C15), we get 



A = 1 + |-^ + 0( ^ ) - (8 ' 26) 



This obviously agrees with Eq. ( 8.25| ) for large N. 

The next quantity we compute is the universal mobility /i* of the QBM model at the intermediate fixed point. 
Firstly, we use the conformal field theory results of the Kondo model. Let us define the spin currents 

J a = v F s z a (8.27) 
= -fd x 4> a = —dt^, (8.28) 

47T 47T 
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and their linear combinations 



J s = 



^ — l^Ja = -7-dt 



j| f = ^O io J a = i-^| f 



(8.29) 
(8.30) 



Note that J s is the total spin current and jf are the spin flavor currents. Since the analogue of the Brownian particle 
position r*j is the total spin flavor in each channel Sf — Oi a S*, the velocity dtTi is mapped to the rate in which 
spin flavor is "injected" into channel i. Since the injection occurs at x = through the exchange coupling Hj, the 
spin flavor injection rate dtSf may be measured via the difference in jf right before and after x — 0. In other words, 
we have a ono-to-one correspondence 



Jf(x = o+) - Jf(x = cr) 



Or 



From the above mapping and Eq. (3.7), we get 

" = ^oi/ dT (1 " e " T )( T - [ J f (°V) - Jf(0-,r)] [jf (0+0) - Jf (0-,0)]) o , 



.31) 



.32) 



where T T is the time ordering operator and the average (• • -)o is calculated with respect to the free non-interactii 
Hamiltonian Hq. The above fixed point correlation functions have been exactly calculated by Ludwig and AffleckBt 
using conformal field theory. Specifically, they are given by 



(T T jf(0 + ,r)jf (0+,0)> = (T T jf(0-,r)jf (0" 0)> Q = 



2t 2 ' 
2r 2: 



where 



(T T jf(0+,r)jf(0-,0)) o = (T T jf(0-,r)jf (0+,0)> = £ 

_ ^1/2/^0 



and the "modular S-matrix" is given by 



iV + 2 



sin 



^1/2/^0 



7r(2j + l)(2j' + l) 



AT • 



Thus, Eq. ( |8.32|) becomes 



1-4' 



= 2 sirf 



iV + 2 



(8.33) 
(8.34) 

(8.35) 
(8.36) 



(8.37) 
(8.38) 



This result agrees with Emery and Kivelsor£3 for N — 2, where the it takes the perfect value /i* = 1 [See also Sec. 
0- 

The perturbation calculation of n* in the large D limit of the QBM model is described in Appendix C. The result 
is given by 



fi* = 2tt — ] + O [ — 7 



N 2 N 3 



1 



AT4 



(8.39) 



It agrees with Eq. (8.3S) for large N. 



21 



IX. ON-RESONANCE CONDUCTANCE AND THE UNIVERSAL SCALING 



In this section, we will compute the on-resonance conductance of the two-lead spin-1/2 quantum dot modelH using 
the mapping of the quantum dot model to both the Kondo problem and the QBM model. For convenience, we will 
actually consider a four-lead spinless model instead of the spin-1/2 model, since they are equivalent (see Sec. [0]). 

First of all, we define the conductance matrix for an TV-lead model 

It characterizes the linear current flowing into the a-th lead in response to an infinitesimal voltage V& in the 6-th lead. 
The mapping described in the previous sections dictates that the current and the voltage are mapped to the spin 
current and the local channel magnetic field in a Kondo model, respectively. Thus, we may alternatively compute the 
"spin conductance" of the Kondo model. Explicitly, it is defined as 

G ^ = -\}^ Q ^TiJ dT(l-e^ T )(T T [j a (0 + ,T)- J o (0- r)] [J 6 (0+ 0)- J 6 (0-,0)]), (9.2) 

where the spin current in each Kondo channel J a is defined in Eq. ( (B.27p . Using Eqs. (B.30) and (8.32), the above 
equation may be evaluated to yield 

Gab = j(jf~ S ab) /*» ( 9 - 3 ) 

where fi is the mobility of the corresponding QBM model. 

A four-lead model is equivalent to the two-lead spin-1/2 model if we map the lead 1 and 2 (3 and 4) to the two spin 
channels of the source (drain) lead. If we apply an infinitesimal voltage V to the source relative to the drain, i.e., 

(Viia = l,2, 

Va -\0 if a = 3, 4, l9 - 4j 

then the current is given by 

4 



I = J2 i-edtQa) (9.5) 

a— 3 
4 4 

a=3 b=l 

= ^V. (9.7) 
Substituting the on-resonance mobility fx* = 1/2 in Eq. ( |S.38 ) to the above equation, we get 

G *=*x = fr < M > 

Physically, this result may be understood as follows: On resonance at T = 0, the conductance of each QPC is given 
by 2/i*e 2 //i = e 2 /h, where there is an overall factor of two due to spin. Then the above G* is the series conductance 
of the two QPCs (see Sec. ||). 

Experimentally, G* may be measured as the peak value of a conductance resonance as one sweeps the gate voltage. 
Again, the shape of the peak follows a scaling function that is independent of microscopic details of the system. Since 
the scaling dimension of the magnetic field Ti is given by 1/3 P the resonance line shape has the form 

G(6no,T) = G(J^Y (9.9) 
In principle, the explicit form of G may be obtained using a Monte Carlo simulation.0 
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X. CONCLUSIONS 



In this paper, we have used a QBM model with periodic potential to study a few different fundamental problems 
in condensed matter physics. Generically, by integrating out unimportant degrees of freedom in a path integral 
representation, one can obtain an alternative description of a problem in terms of a first-quantized quantum mechanical 
model of a single Brownian particle moving in a periodic potential. The remaining degrees of freedom are mapped to 
the spatial coordinates of the Brownian particle, whereas the integrated ones serve as the source of dissipation. 

We have applied this method to such problems as resonant tunneling in quantum dot and Luttinger liquid systems 
as well as multichannel Kondo problems. For an iV-lead quantum dot model in the limit AE <C T <C e 2 /C, we have 
shown that the on-resonance problem is mapped to an (N — l)-dimensional honeycomb lattice and the off-resonance 
problem to an (N — l)-dimcnsional triangular lattice. Each vector component of the Brownian particle position is the 
total charge in each lead, and the periodic potential comes from backscattering. Analogously, by mapping the total 
spin in each channel to each vector component of the position, we could derive similar mappings for the ./V-channel 
Kondo problem in the Toulouse limit. 

The lattice symmetry of the periodic potential determines the zero-temperature phase diagram. For symmorphic 
lattices, there are two trivial zero-temperature phases: If the dissipation is strong or the lattice sites are far apart, 
the Brownian particle gets localized and the mobility /i vanishes. If the dissipation is weak or the lattices sites are 
close to one another, the particle diffuses freely and the mobility takes its perfect maximum value /i = 1. We have 
analyzed the phase diagrams in perturbatively accessible regimes. 

For non-symmorphic lattices, the extra symmetry between sublattices may give rise to an intermediate phase. The 
mobility in this phase is given a universal value in the range < \i < 1 and independent of the microscopic details 
of the problem. The intermediate phase is governed by the multichannel Kondo fixed point. We have explicitly 
demonstrated that the mappings are consistent by computing the scaling dimension of the leading irrelevant operator 
and the fixed point mobility in both the QBM model and the Kondo model, in the large N limit. 

Exploiting the mapping, we have used spin current correlation functions calculated by a conformal field theory to 
compute the mobility at the intermediate fixed point exactly. Using the result 

a* =2 sin 2 - , (10.1) 

we could also calculate the on-resonance linear conductance G a b between each individual lead for the N-\e&d quantum 
dot model. Particularly for the conventional two-lead spin- 1/2 model, the on-resonance conductance is given by 
G = e 2 /2h. Experimentally, this conductance may be measured from the peak value of the conductance resonance 
as a function of the gate voltage. From the QBM theory, the resonance line shape is predicted to be described by a 
scaling function of the form 

G(Sn ,T) = G(^y (10.2) 

where a is determined from the scaling dimension of the local magnetic field at the magnetic impurity. Specifically, 
a = 2/3 for the two-lead spin-1/2 quantum dot model and a = 1/2 for the spinless model. 
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APPENDIX A: RENORMALIZATION GROUP ANALYSIS OF THE QBM MODEL 

In this appendix, we develop a renormalization group analysis method for the QBM model. The renormalization 
group flow equation for the triangular lattice will be explicitly calculated as an example. More examples can be found 
in Appendix C for C-dimensional honeycomb lattices in t he la rge D li mit. 

Let us consider the small v limit of the action as in Eqs. ( |6.l[ ) throu gh (|6.3[ ). It sho uld be straightforward to perform 



a similar analysis for the small t limit using a dual theory in Eqs. (5.10) through (6.12). The partition function is 
given by 
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z = 



J 23[r(r)] 



-So[r]-S„[r] 



^2J-f((-^) n > 



(Al) 
(A2) 



where Zq = J T>[r] e So is the the partition function for v = 0. The brackets denote an average with respect to the 
free action Sq and are defined as 



io) = T Q j " v[v] 0e ~ S °- 



Each term in Eq. (A2) may be explicitly written as 



rr<(-$0 



E 

{G} 



T1<T2<"-<T, 




7rX) 4 G 4 .r(TO 



(A3) 



(A4) 



where vt is a shorthand notation for VG i . The prefactor 1/n! has been absorbed in the explicit time ordering. For 
a non-symmorphic lattice, time ordering is important because the Brownian particle has to move back and forth 
between different sublattices. An explicit calculation at T = shows that 




(Tj _ Ti) 2 +r 2 



if EC = o, 

otherwise. 



(A5) 



Any low energy properties calculated from the above partition function must not depend on the cutoff r c . Therefore, 
by renormalizing r c , we may obtain a set of scaling laws that connect a given problem to another one with renormalized 
paramete rs. In the new theory with cutoff t' c = r c e l , we may define the renormalized amplitud e yf 1 in such a way 
that Eq. (A2) retains the same form when written in terms of vf and t' c . With the help of Eq. (A5), this condition 
may be written as 



vf_ 

T r e l 



(\2 



(roe") 



1 |G ! 



(r 2 -r 1 )2 + (r c e^) 2 



2 r 



(T2-n) 2 +T c 2 



to the second order in v. Assuming r c r c e <C t<l — t%, we get a simple scaling equation 



^edG.r-iy 



(A6) 



(A7) 



Differentiating both sides of Eq. flATf ) with respect to £, we obtain the renormalization group flow equation ( p.9|) . 

Higher order corrections may be calc ulated in a similar way. Let us consider any arbitrary pair of consecutive 
operators at time Tj and r i+1 in Eq. (A4). If r c < T i+1 — < r c e f , they are no longer distinguished in the new theory 
with cutoff r c e l . Instead, they must be treated as a compound operator at a single moment: 



e -t27rG < T(T<) e i27rG i+ iT(T <+1 ) 



i27r(Gi+G i+ i).r(T 4 ) 



(A8) 



This gives rise to the second order correction to Eq. (A7). This process may be thought of as "decimation" of closely 
placed operator sequences in time. In general, we get the n-th order correction by lumping together n operators that 
lie within r c e e . 

In order to show details of the decimation procedure, we will consider a triangular lattice as a specific example. The 
shortest reciprocal lattice vectors Gi and — G, (i = 1, 2, 3) are defined as in Fig. ^(a). We need to find all G-sequences 
that add up to a single reciprocal lattice vector. For example, since Gi + G 2 + G3 = 0, there are two sequences that 
are second order in v: 



{— G 2 ,— G 3 }, {— G 3 ,— G 2 }, 



that add up to Gi. When decimated, they give rise to a correction to the righthand side of Eq. (A7) 

(£{23} + C{32"})^ 2 > 



(A9) 
(A10) 



where 



24 



C{23} - 
C{32} = 



dr' 



dr' 



i2ix -G 2 -r(r 1 )-G 3 -r(r 1 +r')-G 1 -r(T 2 ) 



,i2n -G 3 -r(Ti)-G2T(n+T')-GiT(ra) 



T 2 - Ti 



T2 - n 



2|Gi 



2|Gi 



(All) 
(A12) 



Each of the above terms comes from the corresponding sequence in Eq. ( A9). We have assumed v x = v 2 ~ v 3 = v. 
The overlines in 2 and 3 denote that the corresponding vectors have a minus sign. In general, it is not easy to compute 
C, but if \G t \ 2 w 1, we get 



Thus the scaling equation becomes 



C{23} _ C{ 32 } 



v R e (\Gi\ 2 -m = v + 2iv 2 



(A13) 



(A14) 



By differentiating both sides of Eq. ( A14 ) with respect to I and writing the final equation in terms of th e ren ormalized 
parameter v R , we obtain the renormalization group flow equation. For \Gi\ 2 = 1 + e, it becomes Eq. (7.10). 

This analysis may be straightforwardly generalized to even higher order corrections and nonidentical u's. In general, 
for every sequence {G^, • • • , Gi n } such that Gi p = Gi, there is an n-th order contribution 



C{iil2 ■■■in} V il ' ' ' Vin 



to the righthand side of Eqs. (A7) and (A14), where 

C{ii<2'"i»i} 



0<T2<T 3 <-<T 7l <T c e f \p = 2 



n 



l<p<q<n 



(r q - r p ) 2 + t2 



(A15) 



(A16) 



It is assumed that t\ — in the above integrand. For the two- vector sequences in Eq. (A9), disconnected terms 
identically vanished up to the second order. In general, however, they are finite and have to be carefully taken care 
of. Disconnected terms can be treated basically in the same way as above, except for that they are products of one 
non-loop and one or more loops, where a loop is a subgroup of sequences that add up to zero. For example, there is 
a fourth order disconnected sequence {Gi}{Gi, Ga, G3}. The first subgroup trivially adds up to Gi, but the second 
subgroup is a loop [Fig. 0(b)]. The contribution from this process is given by 



C{l}C{123}Vl«2W 3! 



(A17) 



where we define = 1. Higher order calculations for a D-dimensional honeycomb lattice in the large D limit is 
given in Appendix C. 



APPENDIX B: RENORMALIZATION GROUP ANALYSIS OF THE MULTICHANNEL KONDO 

PROBLEM 



In this appendix, we derive the flow equations ( p. 23 ) and ( p. 24 ) of the multichannel Kondo model. In the boson 
representation, the partition function may be written as 



Z 



V 



WAf]e-! dTH ', 



where H' is the Hamiltonian given in Eq. (8.14). For small Jj_, we may write it in a Taylor series expansion as 



-oE;H#) n ({/fE[^ 

(— TV n (— ) e (e-^M K ^°^) 

n V 2 / Jn<n<-<T n a=l \ T c J M \ 



(Bl) 

(B2) 
(B3) 
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where K = (1 — NJ z /2)/VN and s p alternates between ±1 due to S^ p . The average (• • •) and the free partition 
function Zq are calculated for J± = 0. After the averages are evaluated, the above equation becomes 



N 



where f implies integration over the range n < T2 < • • • < r n , and 



(B4) 



^SpSq 



K 



E°S°* 



(B5) 



It is convenient to view Eq. (B4) as the partition function of a plasma in the r-space. There are N types of unit 
charges O^ 1 (a — 1, • • • , N). Since each charge is characterized by N — 1 vector components (i = 1, • • • , N — 1), it is 
a "vector plasma" . An arbitrary single charge in the plasma may be identified by the position r, the charge type a, 
and the sign of the charge s. Then, V pq is the interaction energy between two changes (r p , a p , s p ) and (r q , a q , s q ). 

We now analyze this model using a real-space renormalization group methodO A similar analysis of ID plasma 
model has been performed for a Luttinger liquid resonant tunneling problem in Ref. [18], but it differs from our 
model in that charges in their model are conventional scalar objects. In the first step of the renormalization group 
transformation, we decimate pairs of closely placed charges. Suppose there is a iO" 1 charge pair at r' and r' + At. 
If At is small compared to the distance between this pair and other charges, the pair may be thought of as a dipole. 
The total interaction of this dipole and all other charges is given by 



K dipol V; At) = -2£ s p s q U 2 + £ 0^0*. At8 t , In ■ 



(B6) 



If r c < At < t c + 5t c , the dipole is no longer recognizable after increasing the cutoff from t c to r c + 5t c . However, 
it induces effective interaction of order 8t c /t c between other charges. When all such dipoles are integrated out, the 
partition function becomes 



n V 7 J a=l \ c / { a? } \ / c p J t p c a 



n „/ N 



a=l {a P } 



■Eh 



™.J 2 z 5t c 2 T q -T f 

' pq r opOqi\ 111 



(B7) 
(B8) 
(B9) 



The notation j means that the sum is only over the vectors that form a loop, i.e., ^2 p s p a ^ — 0. We complete 

the renormalization group transformation by rcscaling r — ► Te e where £ = 5t c /t c . Comparing the final expression 
with the original partition function in Eq. (B4), we obtain the renormalization group flow equations 



dJ± j t n N n 

— - J Z J±{1- -jJ z ), 



(BIO) 
(BH) 



which are exact in J z and perturbative in J±. The extra higher order term in Eq. (8.24) has been read off using the 
mapping between the Kondo model and the QBM model in the limit N 3> 1 and \J Z — 2/N\ -C 1. 



APPENDIX C: PERTURBATIVE CALCULATIONS 



This appendix is devoted to the computation of physical quantities at the intermediate fixed point of the D- 
dimensional honeycomb lattice in the large D limit. In the first subsection, we derive the renormalization group flow 
equation in the Toulouse limit and perturbatively compute the scaling dimension of the leading irrelevant operator. 
In the second subsection, we calculate the fixed point mobility to leading and subleading order in 1/D. 
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1. Renormalization group flow equations 



We use the tunneling part of the action in Eq. (7.281) in the small t limit. The shortest direct lattice vectors of 
a D-dimensional honeycomb lattice are given in Eq. ( |7.38| ). We will continue to use the same convention in which 
Ro G {Ri} ({— Ri}) if h originates from A sublattice (B sublattice). In the Toulouse limit (c = 1), they satisfy 



IRnl 2 



D 



D + l 



(CI) 



If D is large, |Ro| 2 ~ 1 and t is only slightly relevant. The intermediate fixed point thus lies perturbatively close to 
the localized fixed point (t* sa 0). The renormalization group analysis is very simi lar to that of the small v limit in 
Appendix A. The leading order scaling equation may be simply read off from Eq. (A7): 



R„-e/{D+i) 



t n e 



t. 



(C2) 



The next lowest order term is proportional to t 3 and may be computed using a similar renormalization group 
technique as in Appendix A. There are three third-order sequences that add up to R4: 



{Ri, — Ri, Ri}, 
{Ri, — Ri, Ri}, 
{R„ R„Ri}, (i^l). 

According to Eq. ( Al5| ), their contributions to the scaling equation are 

C{iTi}* ' Qiii}^ 3 , ({ill}* 3 , {i l)i 
respectively. There are also disconnected terms 

{Ri}{Ri, — Ri}, 

and their contributions are 



{«}• 



Using Eq. (|A16| ) , we may explicitly compute the sum of all third order corrections: 



D+l 



Cs - C{ni} + ^2 + _ ^2 ^{^{ii} 



D+l 



2t / dx 



D 



{\-x) 2 ' D -1 



In the large D limit, 



Cs = -(D + 5)£ + 0{l/D). 



(C3) 
(C4) 
(C5) 
(C6) 

(C7) 
(C8) 

(C9) 



In order to obtain results up to the subleading order in 1/D, we need to calculate the fifth order contribution £5. 
A calculation analogous to the above yields 

Cs = C{ 11111} 

+ X] (QiTiii} + Quia} + QiiiTi} + QiiiTi} + QiTlii} + QiiiTi] 

i>2 



\C{Uiii} ^{iilii} + ^{iiiil} 



2<i<j 
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C{ 1 1 1 } + 2J (Qim} + Qiil} 



i>2 



C{ii} 



i>l 



Cm 



Cm 



X! Qiiii} + X! (Qiijj} + C{yji} 
«>1 l<i<j 
2 



C {J 7 } 



Collecting only the leading order terms in D, we get 

C 5 - 2D 2 £ + 0{D). 



(CIO) 



(Cll) 



Combining Eqs. (C2), (|C9|), and (Cll) and differentiating it with respect to £, we obtain the following renormalization 
group flow equation: 



^ = -L-t R -(D + 5)(t R f + 2D 2 (t«f. 



From dt/di — 0, the intermediate fixed point is given by 



By substituting t R = t* + St in Eq. flC12| ), we get 

dSt 2 6 



0( 



1 



dt D D 2 y D 3 '' 
Therefore, the scaling dimension of the leading irrelevant operator is 



(C12) 



(C13) 



(C14) 



(C15) 



Since the Z?-dimensional lattice model is mapped onto the (D + l)-channel Kondo model, we may compare Eq. (C15) 
with Eq. ( p.25| ) by setting N = D + 1. Clearly, they agree with each other. 



2. Fixed point mobility 

Now we compute the mobility at the intermediate fixed point obtained above. We will first add a source term 

1 



S h = I / duj |w|e MT <= [h(w) • k(-w) + c.c] 



to the total action. From Eq. (|6.13 ), the mobility may be expressed as 

1 ° f 

fi=l— lim — : — r } / dui' ujlu' Cki(u))ki(—uj')) 
<»-+oD\(j\^J 

1 ^ 1 f , / / s 2 z 

= 1 - hm — — r > — du> LULU — — — , 

h,a>->o D\lu\ ^ Z J 5hi(-uj)Shi(u)') 



where 



Z = J £>[k(r)] e- s[k] 
is the partition function. By shifting the dummy variable 



(C16) 

(C17) 
(C18) 

(C19) 
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k(cj) -> k(w) - h(w), 

we may write the action as 

S= l -Jdu MeM T « [|kH| 2 + |hH| 2 ] £ 



3 -i27rRi-(k-h) 



H.c. 



Now it is easy to evaluate Eq. ( |C18| ): 

2tt 



/i = lira 



E 



R R' /" dr 



R,R' 



— 1-e" 
7V 



^2 / i27r[R-k(r)+R'-k(0)]\ 



(C20) 



(C21) 



(C22) 



where (• • •) denotes the connected part of an average. In the small t limit, we may Taylor series-expand the above 
equation and write 



/ - +"+2 , 
t 2 / e i2-7r[R-k(r)+R'-k(0)]\ _ ^ l _ / e »2»r[H..k(T)+R/-k(0)] (-g t ) 



m+2 



i=l 



3 i27r[R-k(T)+R'-k(0)+^. Ri-kfa)] 



(C23) 
(C24) 



where (• • -)q is the connected part of an average over So [Eq. ( 6.1 1| ) ] . Note that the lattice vector is alternately chosen 
from {Ri} and {— R;}. The first non- vanishing term is second order in t and is given by 



M (2) = Hm K / drfl - ) / e ^R-[k(r)-k(o)]\ 

w^o \u\t% ^ D J v '\ /o 

Ft 



2nt 2 



= lim r — - / dr(l - e 4 

= 27r 2 t 2 (|^|r c )- 2 /( Z3+1 ). 



From a similar calculation, the next order term is given by 



D/(D+l) 



M W = Hm JE* y / ^^^(X.^rW^lRkM+R'-kTO+R.-kCrO-fR+R'+RO-k^)]^ 

^"M^^ D J \ /o 



lim 



4vrL>t 



4 /■/ 



drdri dr 2 



1 



1^1 •/ " "[(r 2 -T 1 ) 2 +T c 2 

= -47r 2 £>t 4 [l-ln(Mr c )], 

where we kept only the highest order terms in D. In Eq. ( p29| ), J means an integration over the range 

ri < < r 2 < r and < ti < r < t 2 . 



(C25) 

(C26) 
(C27) 

(C28) 

(C29) 
(C30) 

(C31) 



Combining Eqs. ( |C27 ) and (C30) and writing them in terms of the renormalized tunneling amplitude t R , the mobility 
is given by 



fi = 2tt 2 



{t R f-2D (t 1 



Substituting the fixed point value t* in Eq. ( C13| ), the above equation becomes 



D 4 



(C32) 



(C33) 



When N = D + 1, this is consistent with the exact solution Eq. (p.38). 
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D 



lattice {R} 



|Go| : 



dv/dl 



triangular 
honeycomb 



l + e 
1 - e 



-ev + 2v 2 
ev - 2v 2 



1 - e 
l + e 



3jt± 2 

3tt 2 2 
2 fc 



FCC 
diamond 



l + e 
1 - e 



-ev + (A 3 + B 3 )« 3 y 
ev - (A 3 - B 3 )v :i £ 



1 - 2e 
1 + 2e 



16ti-^ 
3(A 3 +B 3 ) 

16tv 2 
3(A 3 -B 3 ) 



D > 4 



generalized FCC 
generalized diamond 



l + e 
l + e 



-ev + B D v 6 
-ev + B D ?) 3 



1 -2e 
1 -2e 



1 _ 4(D + 1)^ 

DBfl fc 
1 _ 4(D + 1)^ 2 



^ 3 = 6B(l/3, 1/3) = 6 r(1 ^W «l « 31.80 

B D = 2£ £ P^T^T- 2 + (t^) 2/D - ll , (Bs = 4737T « 21.77, B x = 0) 



TABLE I. Results of the renormalization group analysis near |Go| 2 = 1 for small v. We assume e > 0. fx* is the mobility 
and A is the scaling dimension of the most relevant operator, both at the fixed point v*. 
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3^2 
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et - 3t 3 
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1 + 2e 


n 2 e 
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FCC 




l + e 


-et + tt 2 


e 
4 


1 - e 


& 




diamond 




1-e 


et - C 3 t 3 


\FSi 


1 + 2e 


3C 3 £ 




generalized FCC 
generalized diamond 


l + e 
1-e 


-et + 2(D - l)t 2 
et - C D t 3 




1 - e 
1 + 2e 


(D+lV 2 


D > 4 


2(D-1) 

1 € 

V c D 


2(£>-l) a 6 

2(D + 1)7T 2 

DC D 


Co 


= D + l-2f t 'dx^ 


D(l- 


x f/o — D + 2x] 


, (C'2 = 3, Coo — D + I 


+ 0(D~ 1 )) 







TABLE II. Results of the renormalization group analysis near |Ro| = 1 for small t. We assume e > 0. fx* is the mobility 
and A is the scaling dimension of the most relevant operator, both at the fixed point t* . 
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2DEG 




FIG. 1. Schematic view of a quantum point contact formed in a two-dimensional electron gas (2DEG). The shaded regions 
are metal gates and the thin lines are the boundaries of the 2DEG. Negative voltage on the gates depletes electrons underneath 
them and separates the 2DEG into two regions connected only by a narrow channel. 




FIG. 2. Analogy between a quantum point contact and the quantum Brownian motion in a tilted periodic potential. 




FIG. 3. Minima of the periodic potential in Q1-Q2 space. After n has been integrated out, the on-resonance model is 
described by the corrugated ID lattice. The two sets of degenerate sites in the lattice are depicted by filled circles and squares. 



(a) (b) 




FIG. 4. (a) Schematic view of a three-lead quantum dot model, (b) Two adjacent planes of triangular lattices, or equivalently, 
a corrugated honeycomb lattice, which are embedded in a cubic lattice. It is a view from the (111) direction. 
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(a) 



(b) 
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—I -1 ^ 



1/2 
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1 IR/ 



2D/(D+1) IG I 

T 1 ► -1 V 



+ 



4- 







ir i 



2D/(D+1) 1 

FIG. 5. Flow diagrams for symmorphic lattices: (a) triangular lattice; (b) D-dimcnsional generalized FCC lattice. The top 
(bottom) line in each diagram represents the small v (t) limit. Stable (unstable) fixed points are depicted by solid (dotted) 
lines, and arrows indicate the renormalization group flows. The fixed point mobility is known exactly at the self-dual point 
(dashed line) and perturbatively near |G„| 2 = 1 and jR j 2 = 1. 



O 



O 



-R, 



o 



-R, 



O 



o 



A sublattice 
O B sublattice 



• • • 

FIG. 6. Shortest displacement vectors in a honeycomb lattice. 




(ve'f) (ve-'f ) 3 = vV 1 * 

FIG. 7. (a) Reciprocal lattice vectors of a honeycomb lattice are shown with the amplitudes, (b) A phase factor e ±ln is 
acquired from a trip around a triangle, which is analogous to a tight-binding model with flux ±7r per plaquette. 
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4/9 |R Q | 2 



(b) 




i \G r 



i v- 



1/2 



9/16 |R Q | 2 



(D/CD+l)) 2 D/(D+1) 




(c) 



t 

1 D/(D+1) (D/(D+1)) 2 |r/ 
FIG. 8. Flow diagrams similar to Fig. 4 for non-symmorphic lattices - (a) honeycomb lattice, (b) diamond lattice, and (c) 
their generalizations to D > 4. The fixed point mobility is known exactly at the Toulouse limit (dashed lines) and perturbatively 



near |G 



o 



1 and |R, 




FIG. 9. Schematic diagram of renormalization group flows for the distorted honeycomb lattice model, which describes 
resonant tunneling of spin- 1/2 Luttinger liquid [Ref. 0. 





2/N 
FIG. 10. Flow diagram of the N channel Kondo model for small Jj 
strong coupling fixed point is marked with the full circle. 



4/N J z 

The dashed line is the Toulouse limit, J z = 2/N. The 
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